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Abstract. In this article we improve some of the inverse spectral results proved by Guillemin and Uribe 
in |GUI . They proved that under some symmetry assumptions on the potential V(x), the Taylor expansion 
of V{x) near a non-degenerate global minimum can be recovered from the knowledge of the low-lying 
eigenvalues of the associated Schrodinger operator in R™. We prove some similar inverse spectral results 
using fewer symmetry assumptions. We also show that in dimension 1, no symmetry assumption is needed 
to recover the Taylor coefficients of V(x). We establish our results by finding some explicit formulas for 
wave invariants at the bottom of the well. 

Keywords: Schrodinger operator, Wave invariants, Semi-Classical trace formulaes. 

1. Introduction and Statement of Results 

In this article we study some inverse spectral problems of the eigenvalue problem for the semi-classical 
Schrodinger operator, 

h 2 

(1) P = ~A + V{x) on L 2 {R n ), 
associated to the Hamiltonian 

Here the potential V(x) in fl} satisfies 
( V(x) G C°°(R n ), 

(2) < V(x) has a unique non-degenerate global minimum at x = O and V(0) = O, 
[ For some e > O, V _1 [0, e] is compact. 

Under these conditions for sufficiently small H, say H G (0,ho), and sufficiently small S, a classical faet tells 
us the Spectrum of P in the energy interval [O, S] is finit e. We denote these eigenvalues by 

{E 3 {h)YjU- 

We call these eigenvalues the low-lying eigenvalues of P. We notice the Weyl's law reads 

(3) rn = N H (å) = tt3lO<E j (h)<6} = -±—([ dxd£ + o(l)). 



Recently in |GUj , Guillemin and Uribe raised the question whether we can recover the Taylor coefficients 
of V at x — O from the low-lying eigenvalues Ej(h). They also established that if we assume some sym- 
metry conditions on V, namely V(x) = f{%\, then the I-parameter family of low-lying eigenvalues, 
{Ej(H) | h G (O, ^o)}> determincs the Taylor coefficients of V at x = 0. 

In this article we will attempt to recover as much of V as possible from the family Ej(H), by establishing 
some new formulas for the wave invariants at the bottom of the potential (Theorcm ll.l[) . Using these new 
expressions for the wave invariants, in Thcorcm 11.21 we improve the inverse spectral results of |GU| for a 
larger class of potentials. 
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A classical approach in studying this problem is to examine the asymptotic behavior as h — > of the 
truncated trace 

(4) Tr{e(P)e^ Lp ), 

where O G Cq°([0, oo)) is supported in / = [0, S] and equals one in a neighborhood of 0. 

The asymptotic behavior of the truncated trace Tr(t,h) around the equilibrium point (x,£) = (0,0) has 
been extensively studied in the literature. It is known that (see for example |BPUj ) for t in a sufhciently 
small interval (0, to), Tr(0(P)e^~ p ) has an asymptotic expansion of the following form: 

oo 

(5) Tr(Q(P)e = i Lp )^^2a j (t)h j , h^0. 

j=o 

Throughout this paper when we refer to wave invariants at the bottom of the well, we mean the cocfhcients 
a 3 (t) in ©. 

By applying an orthogonal change of variable, we can assume that V is of the form 

1 - 

(6) V(x) = -J2"l4 + W(x), u; k >0, 

L k=i 

W(x) = 0(|x| 3 ), |x|->0. 

In addition to conditions in ©, we also assume that {w/.} are linearly independent over Q. We note that 
we have W(0) = VVK(0) = HessW(O) = 0. 

Our first result finds explicit formulas for the wave invariants. 
Theorem 1.1. For < t < mini< fe <„{ 2^} ; 

n 1 li" 
a Q (t)=Tr(e^r H °) = ]J — — , where H = --h 2 A + - ]T 

k=i Zlsm ~ 1 1 fc=i 

2. For j > 1 , the wave invariants aj(t) defined in |5|) are given by 

(7) aj{t) = aorøy" i ' ( "" 1)+f ellS3nHl f 1^ - C " Pi+M0)dsi...ds u 

l=1 Jo Jo Jo 

where for every m, 

i—m 

PnMO) = < Hf l V,V > m (6,)(0), 

2—1 

and i?T is i/ie inverse matrix of the Hessian Hi =Hess ^>i(0), where 

n l 
k=l z=l 

The Hessian of is calculated with respect to every variable except t. Therefore the entries of the matrix 
Hf 1 are functions of t. The matrix Hf 1 is shown in 



1. 
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3. The wave invariant o,- (i) is a polynomial of degree 2j of the Taylor coefficients ofV. The highest order of 
derivatives appearing in aj(t) are of order 2j + 2. In faet these higher order derivatives appear in the linear 
term of the polynomial and 
(8) 

a j(t) = (2°V+ 1 ~( ~2 ^ a ^25 +2 ^(^) + i 0, polynomial of Taylor coefficients of order < 2j + 1} 



Notice that in |JS]), we have used the standard shorthand notations for multi-indices, i.e. d — (arj., ...a n ), 



tf = (cJi, |a| = a\ + ... + a n , a\ = a\\...a n \, X a = X" 1 —X% n , and D™ = ni r Q „ with m = 



Our second result improves the result of Guillcmin and Uribe in |GUj . This theorem is actually a non- 
trivial corollary of Theorem 11.11 

Theorem 1.2. Let V satisfy (0, (0), and be of the form 

(9) V(x) = f(xl...,x 2 n )+x 3 n g(xl...,x 2 n ), 

for some f,g £ C°°(R n ). Then the low-lying eigenvalues of P = — ^h 2 A + V determine D^V(O), \et\ — 2,3, 
and if D^g n V(0) := §^r(0) ^ 0, they determine all the Taylor coefficients ofV at x = 0. 



One quick consequence of Theorem 11 .21 is the following: 

Corollary 1.3. Ifn — 1, and V € C°°(K) satisfi.es ([£[). then (with no symmetry assumptions) the low-lying 
eigenvalues determine V"(0) and T^( 3 )(0), and ifV^ 3 \0) ^ 0, then these eigenvalues determine all the Taylor 
coefficients of V at x — 0. 

Let us briefly sketch our main ideas for the proofs. First, because of a technical reason which arises in 
the proofs, we will need to replace the Hamiltonian P by the following Hamiltonian H 

( H(x,t) = %e + V h (x) i 

Wfi(x) = x( i )W(x), e > sufficiently small, 

where the cut off \ e Co°(R n ) is supported in the unit ball -Bi(O) and equals one in Bi (0). 

Then in two lemmas (Lemma 12. II and Lemma [2~2]) we show that for t is a sufficiently small interval (0, to), 
in the sense of tempered distributions we have 

Tr(&{P)e^r p ) = Tr(e^ É ) + 0(h°°). 

This reduces the problem to studying the asymptotic of Tr(e~ H ). For this we use the construction of the 
kernel k(t,x,y) of the propagator U(t) = e~ H found in (Zj. In [Z] it is shown that 



(10) k(t, x, y) = C(t)ei s ^ ]T ai(t, ti, x, y), 

1=0 

where 

" uj 1 

(11) S(t,x,y)=J2^ jL -A~(cosu k t)(xl+yl)-x k y k ), 



k—l 



ao = 1, and for l > 1, 



a l (t,h,x,y) = (^) ln (±) l ( n+ V C ... f ... f e^%{s,x, y ,z,i)d l zd^d [ s, 
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where 

n / 

®i = E{ y cot w **(^*) a + £(4-i - 

k=l i=l 

and 

i. -TT™- f COSU} k Si ^ k k waMkSj k smuj k (t - s t ) n k sina^Sj fc 

°i I I o \ z i+l + z i ) Si i : 7 V ~l : 7 1 J- 

2 — 1 



Next we apply the expression in (fT0|) for fc(t, x, y) to the formula Tr(e~ H ) — J k(t, x, x)dx. Then we ob- 
tain an infinite series of oscillatory integrals, each one corresponding to one etj. Finally we apply the method 
of stationary phase to each oscillatory integral and we show that the resulting series is a valid asymptotic 
expansion. From the resulting asymptotic expansion we obtain the formulas 0. 

Now let us compare our approach for the construction of k(t,x,y) with the classical approach. In the 
classical approach (see for instance |DSj| ,[D]. [R] . |BPU] and [U]), one constructs a WKB approximation for 
the kernel kp(t,x,y) of the operator <d(P)e~ p , i.e. 

(12) k P (t, x,y) = J e^(t^n)-y-n) bp ^ Xj y ^ ^ h)djh 

where (pp(t,x,r)) satisfies the Hamilton-Jacobi equation (or eikonal equation in geometricai optics) 

d t (pp(t,x,rj) + P(x,d x (p P (t,x,r))) = 0, <pp\t=o = x.rj, 
and the function bp has an asymptotic expansion of the form 

oo 

b P (t,x,y,ri,h) ~ ^2b Pij (t,x,y,Ti)h j . 

3=0 

The functions bp t j(t, x, y, 77) are calculated from the so called transport equations. See for example [R], |DSj| , 
[EZj or Appendix A of the paper in hand for the details of the above construction. 

In this setting, when one integrates the kernel k(t,x,y) on the diagonal and applies the stationary phase 
to the given oscillatory integral, one obtains very complicated expressions for the wave invariants. Of course 
the classical calculations above show the existence of asymptotic formulas of the form {SJ (which can be used 
to get Weyl-type estimates for the counting functions of the eigenvalues, see for example [BPUj V Unfortu- 
nately these formulas for the wave invariants are not helpful when trying to establish some inverse spectral 
results. 



Hence, one should look for more efhcient methods to calculate the wave invariants aj(t). One approach is 
to use the semi-classical Birkhoff normal forms, which was used in |GUj . to get some inverse spectral results 
as mentioned in the beginning of the introduction. The Birkhoff normal forms methods were also used before 
by S. Zelditch in [Z4] to obtain positive inverse spectral results for real analytic domains with symmetries 
of an ellipse. Zelditch proved that for a real analytic plane domain with symmetries of an ellipse, the wave 
invariants at a bouncing ball orbit, which is preserved by the symmetries, determine the real analytic domain 
under isometries of the domain. 

Recently in |Z3| , Zelditch improved his earlier result to real analytic domains with only one mirror symme- 
try. His approach for this new result was different. He used a direct approach (Balian-Bloch trace formula) 
which involves Feynman-diagrammatic calculations of the stationary phase method to obtain a more explicit 
formula for the wave invariants at the bouncing ball orbit. 

Motivated by the work of Zelditch [Z3j mentioned above, our approach in this article is also somehow 
direct and involves combinatorial calculations of the stationary phase. 

Our formula in (JTUJ) for the kernel of the propagator, U(t) = e~ H , is different from the WKB-expression 
in the sense that we only keep the quadratic part of the phase function, namely the phase function S(t, x, y) 
in (jlip of the propagator of Anisotropic oscillator , and we put the rest in the amplitude ^2^QO,i(t,h,x,y) 
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in ITU1). The details of this construction are mentioned in Section 2.2. 



Remark 1.4. After the initial posting of this article, Guillemin and Colin de Verdiére posted two articles 
(see [CGlj . also [C],) in which they study some inverse spectral problems of 1 dimensional semi-classical 
Schrddinger operators. One of the main results in |CGlj is our Corollarv \1.3\ in this paper. 

1.1. Acknowledgements: I am sincerely grateful to Steve Zelditch for introducing the problem and many 
hclpful discussions and suggestions on the subject. I would also like to thank him for his great support and 
encouragement as I was writing this article. 

2. Proofs of the results 

2.1. Some reductions. Because of some technical issues arising in the proof of Theorem ll.li we will need 
to use the following lemmas as reductions. 

In the following, we let \ € C^°(R") be a cut off which is supported in the unit ball -Bi(O) and equals one 
in 51(0). 

Lemma 2.1. Let the Hamiltonians P and H be defined by 

( H{x,£) = \e + V h {x), 

(13) 



Wn,(x) = x( i_ )W(x), £ > sufficiently small, 

h2 e 



and let P and H be the corresponding Weyl (or standard) quantizations. Then for t in a sufficiently small 
interval (0, tg) 

Tr{Q{P)e^r p ) = Tr(Q(H)e^ å ) + 0{h°°). 
In other words, the wave invariants aj(t) will not change if we replace P by H. 
PROOF. Proof is given in Appendix A. 

Next we use the following lemma to get rid of <d(H). 
Lemma 2.2. Let H be defined by ((73J). Then in the sense of tempered distributions 

Tr(Q(H)e^ 1É ) = Tr(e^) + 0(R°°). 
This means that if we sort the Spectrum of H as 

Ei{K) < E 2 (h) < ... < Ej{K) -» +oo, 
then for every Schwartz function tp(t) G S(M) 

< Tr(e^ É ) - Tr(e(H)e^ É ), <p(t) >= £(1 - Øførø))0(M&) = 0(ft~). 

PROOF. Proof is given in Appendix B. 

Because of the above lemmas, it is enough to study the asymptotic of Tr(e^~ H ). 
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2.2. Construction of k(t,x,y), the kernel of e » H . In this section we follow the construction in [Z] to 
obtain an oscillatory integral representation of k(t,x,y), the kernel of the propagator e~ H . The rcadcr 
should consult Z] for many details. In that article Zelditch uses the Dyson's Expansion of propagator to 
study the singularities of the kernel k{t,x,y). But he does not consider the semi-classical setting h — > in 
his calculations (i.e. h = 1). So we follow the same calculations but also consider h carefully. 
In [Z] , potentials of the form 



1 ™ 

k=l 

are considered, where W% S B(M. n ), i.e. bounded with bounded derivatives. 
We denote 

H = — ifi 2 A + | J2k=i ^fc^fc' (Anisotropic Oscillator) 

H = H + W H (x) = -ift 2 A + V H {x), 
and by Uo(t) = e~ H °, and U(t) = e~ H , we mean their corresponding propagators. 

From 

(iHd t -H )U(t) = W h .U(t), 

we obtain 

1 '* 



(14) U(t) = U {t) + - / U (t - s).W h .U{s)ds. 

™ Jo 

By iteration we get the norm convergent Dyson Expansion: 



(15) U(t) = U (t) + V ttL f / ' 1 U (t)[U (s 1 )-\W H Msi)}AUo(si)-\W H .U (s l )}ds l ...ds 1 . 
1=1 \ ' Jo J° 



For t ^ — , the kernel of U (t) is given by 



(16) k Q (t, Xl y) = (]]_— -± -)i e *S(t,*,v) 

J - J - 2irinsmujkt 

k=l 



U) 1 

S(t,x,y) = V -(-(cost^)^ +?v|) - x fc y fe ). 



where 

n 

S'(/ X. 7/1 = V 

sinwfct 2 

1 

Then by taking kernels in (fT5|) and after some change of variables (see [Z]), we get 

n oo 

(17) k(t,x,y) = (TJ — -)M s ^5>(i,M,2/), 

fe=i ;=o 

where ao = 1 and for i > 1, 

2/ 



(18) oi(*,».^») = (2^)' n (^)' (n+1) / - / / - I <**'b,[.s..r. !l .z:.C)<l , z<l , l,<l'». 

where 



o Jo 



si-l 



n l 

(19) *' = E< Y cot Wfct(zife)2 + " 

fc=l i=l 
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and 



20 bi=\\Wh\ {Zi+i+Zi) Ci + -■ 7 V +—. rx ). (z i+ i:=0) 

z— i 

In [Zj, integration by parts in the variables z, £, for the integrals in (|18| are performed (in that article, there 
is no h, i.e. K= 1, because of different motivations) to prove that the sum a(t, 1, x, y) = X]j=o a z(^; 1> æ ! 2/) ^ s 
absolutely uniformly convergent. Additionally, it is shown that a(i, 1, as, y) is in 5(M™ x K"). More precisely, 
there exists ko = ko{a, (3,n) such that 

(21) |^ 0l (t,l,x,tf)| < ic QAn (t) Z ||W 1 ||{ Q|+| ^ |+feo , (Wi = W h \ h=l ) 

which implies 

|c>«c?f«(i, l,cc,y)| < e X p{C Q ^, ra (t)||W 1 ||| Q | + | /3 | +fco }- 
The estimates (|21l) will change if one considers h in the calculations. We will establish these h estimates in 
Lemma |2~51 As a simple consequence of Lemma 12. 31 let us assume for now that a(t,h,x,y) € -B(R™ x K"). 

2.3. Oscillatory Integral Representation for the Trace of U(t) — e^^ 11 . In this section we show that 
the integral TrU{t) — J k(t, x, x)dx is convergent as an oscillatory integral and that from (fl"7f and (p~8|) we 
can write 

(22) 

21+1 

TrU(t) = (f[ ^ )lf](-l)'"(l)'("+i) /"*... / e^^^x.x.^eV^^d'^. 

2mhsmuj k t 2rr iTi Jo Jo J J 

whcre 

n / 

(23) * z = x, *) + $,= £{(-w fc tan 4 + ^ cot w *t(«f) a + $>*_i - zf )£*}. 

fe=i j=i 
Before proving (f2"2"| . we review some standard facts. First of all we know that the sum 

TrU{t) = 2je i— 

is convergent in the sense of tempered distributions, i.e. TrU(t) G S'(M). This can be shown by the Weyl's 
law in its high energy setting, which implies that for potentials of the form V(x) = | Sfe=i ^fc^fc + Wft(x), 
with W € B(M. n ), for fixed H, the j— th eigenvalue Ej(H) satisfies 

(24) Ej(H)~C(n,H)j*, j^oo. 
Another way to define Tr £/ (i) is to write it as the limit 

e R . 

e^w o-ru. 

This time the Weyl's law ([24| implies that the sum Tr U(t — ie) is absolutely uniformly convergent because of 

the rapidly decaying factor e s . As a result, U(t — ie) is a trace class operator. It is clear that the kernel 
of U(t — ie) is k(t — is, x, y), the analytic continuation of the kernel k(t, x, y) of U(t). Clearly k(t — ie, x, y) 
is continuous in x and y. So we can write Tr U(t — is) = J k(t — ie, x, x)dx. We notice that this integral is 
uniformly convergent. This is because up to a constant this integral equals to J e'E 5 ^' -46 '") a(t — ie, h, x, x), 
and the exponential factor in the integral is rapidly decaying for e > as |x| — > oo and a is a boundcd 
function. More precisely 

RØS(f - ie, x, x)) = ± tan(^p£l))4 = £ "f~^ 2 xt, 

fc=i fe=i 1 1 
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and 



< 0. 



|1 -f e w k (it+e)\2 

The discussion above shows that the integral J k(t, x, x)dx can be defined by integrations by parts as follows: 
Sincc 



< D x > 2 e iS ^ x ^ := (1 - A)e^'^ = (1+ || 2£tan(— )x || 2 +2^ £ cu k tan(^))e lSl ' ' • ' ! 

k=l 



we can write 

(26) / e* s{Lx < x) a(t,h,x,x)dx = ti* I e iS{t > x < x) a(t,h,Vhx,Vhx)dx = 



R f /" e iS(t,x,x) ^ >2 2c Jtan(— )x || 2 +2^u fe tan( — ))~ 1 )"°a(t, Vftx, Vfcc)<£c 

^ ~ fc=i 

If we assume < t < mini<fc<„{ .^j-}, then by choosing tiq > and because a(t, H, x, y) G 2?(R™ X R"), the 
integral becomes absolutely convergent. Finally, since by (|27j) the series a(t,H,x,y) = ~}2dL Q ai{t,h~,x : y) is 
absolutely uniformly convergent, we have 

/ e^&^otøM. = / e*^*'* 1 ^ oi(t, ft, se, x)dx, 
•* i=o J 

and therefore we obtain (|22p . which is an infinite sum of oscillatory integrals. The next step is to apply 
the stationary phase method to each integral in (|22|1 and then add the asymptotics to obtain an asymptotic 
expansion for the TrU(t). Because we have an infinite sum of asymptotic expansions, we have to establish 
that the resulting asymptotic for the trace is a valid approximation. Hence we have to find some appropriate 
Æ-cstimatcs for the rcmaindcr term of the series 



2.4. ft-estimates for the remainder. The goal of this section is to find an ft-estimate for the remainder 
term of the series (|22[) . First we state the following important estimate, which shows how the estimates in 
(f2"Tj) change in the ft-dependence case. 

Lemma 2.3. For every a, (3, there exists ko = ko(a, (3) such that for every < h < ho < 1 

\ l \lUf*\\l 



(27) \d^ ai (t,H,x,y)\ < C ^»®^\\ W+ W+* M-*. )t 

where 

. . i i . Wnifrhx) ... .W(hix) 

is uniformly in B(W™); i.e. Wj* is bounded with bounded derivatives and the bounds are independent of h. 

Before proving this lemma, let us show how the lemma is applied to get estimates for the remainder term 
of the series (|22j) . Define Ii(t,h) to be the i-th term in (|22| (removing the constant Y[k=i( 2Trisinui k t )^ ) ' 

(29) I I (t,R) = fi-9 / e ^ s{t ' x ' x) ai(t,h,x,x)dx. 



Hence by this notation, TrU(t) = J\k=i ( 2^ sinoj fc t )^ S/^o^(^ fy- If m p9|) we integrate by parts as we did 
in |26|) . and choose n = [§] + 1, then using (j2"7)) we get 

m,n)\ < g;(t) c ' w( * ),|l ^" r|l ^ +fc ° fi f (^- ae ), 

where C n (t) — max\ a \+\p\<2n {Ca,p,n(t)}- We choose e > such that i — 3e > 0, or e < \. Now it is clear 
that for every positive integer to, and every < h~ < ho < 1, 
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(30) \"jTli(t,h)\ <C' n {t)e^ c ^ w ^o+^h m ^-^. 



Since by Lemma \2 .3[ sup 0<R<1 ||Wjj||2no+*io < 00 > we have 



(31) TrUit) = JJ(_^_)é £ ^(t.fi) + 0(ft«^-^). 

fc=l fc ; =0 

This is a useful estimate that we will use in the next section to obtain the asymptotic expansion of the trace 
and prove Theorem ll.il Let us first prove Lemma [ 



Proof of Lemma [2731 

The proof is straightforward from (|2ip . First, in (|18p . we apply the change of variables z i— > h^z and 
C i — * ^2^. This gives us h ln in front of the integral. Then we replace W% by h 3 ^ e ^W£. After collecting all 
the powers of h in front of the integral we obtain 

21 

<H(t, h, x, y) = {ZL) ln h 1 ^-^ j ... J ... | é*%{s, x, y, z, £)d l zd l ^d l s, 

where 

f-^ 2 Wfc smut! smwfci 

z— 1 

Next we apply (f2"Tj) to the above integral with W\ replaced by W£, and we get (|2T|) . To finish the proof we 
have to show that for every positive integer m we can find uniform bounds (i. e. independent of K) for the 
m-th derivatives of the function W£(x). Since x( x ) is supported in the unit ball, from the definition (|28|) we 

see that W£ is supported in \x\ < h~ e . So from ([28]) it is enough to find uniform bounds in h for the m-th 

i 

derivatives of the function W } h \ 2x ] in the ball la; I < h~ e . This is verv clear for m > 3. For m < 3 , we use 

H Hi-e) ii ■' — 

the order of vanishing of W(x) at x — 0. Since W(Q) — VT / F(0) =HessVF(0) = 0, the order of vanishing of 
W at x = is 3. Therefore in the ball |x| < h~ £ , the functions 

W(h^x) d a W(h?x) d a d^W{fåx) 
{hhx) 3 ' (h^x) 2 hhx 

are bounded functions with uniform bounds in h, and the statement follows easily for m < 3, noting that 

\x\ < fr £ . 

2.5. Stationary phase calculations. In this section we will apply the stationary phase method to each 
Il(t,H) in d3l|). We know 



2/+1 



Si-l 



(32) Ifah) = (-i)'"(4) K " +1)+f / .» / / - / e^hisMZ^Ssd'zd^dx. 



2tt ih 

It is easy to see that the only critical point of the phase function VE^, given by (f2"3")l , is at (x, z,£) = 0. 
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Next we calculate 77; =Hess>I';(0) and 77 ; 1 . In the following we use the notation D(v) for the diagonal 
matrix Diag{v\, v n ), where v = {v\, ...,v n ). From (f2"3")l . we get 



(33) 


/ £>(-2c3tan(f ))„ x „ 
























( -i o 







/ D(ujcotujt) nxn \ 

l o o) 




7-/0 
I ... 






In X In 




Hi = 




f-I I .. \ 
-I I .. 




V o 







■ o ... ... 

/ 









V 


V o o -i j 


In x Hn 




whcrc / = I nxn is the identity matrix of size n x n. 










' K \ 






Since Hi is of the form Hi = 


A B , the inverse matrix cquals iZ, 1 






v B T 







o \ 
o 



1 1 1 



\ 



In x In 





) 


(2/4 


l)nx(2i+l) 


K- 1 


















b t1 





B- 1 - 


-B- X AB T 



A simple calculation shows that 

/ ^cot(f)) 



(34) 



o 



V 

where f2 = D(w cotujt). 
It is also easy to see that 



/ -I ... 

-I -I 

-I . 

\-I ... 



-/ o 

1 - 1 ) 










( 


-J -1 ... 


-/ \ 




-I -1 ... 


-/ 




0-/ ... 








-/ 


\ 


... 


1 ) 



n — ri 



-n 
-n 



\ -o -n ... -n / 



(35) 



det/7; = (-!)('+!)" Y[ 2w fc tan 



fc=i 



UJkt 

2 



By applying the stationary phase lemma to (132| and plugging into (|31j) we obtain 
(36) 



n m—l — l 



(— iy n e i2 i sgnHl 



,-/(n+l)+4 



v/|det// ; | 



o Jo 



^ tv> P j b l (0)ds l ...ds 1 + 0{h m ^-^)) 

3=0 



where 



(37) PM*,W 



' J ' r lr ..,r 2j 6A i "' ZJ 



where in the sum (|37p the indices 7*i,...,r2j run in the set .A; = { 
r, r' (z Ai, corresponds to the (r, r')-th entry of the inverse Hessian HT 



X . Z\ . ..z 



and fcp' with 
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We note that Pj-6j(Q) = if 2j < 31. This is true because of {H} and because W(0) = VW(0) = 
HessW^O) = 0. This implies, first, there are not any negative powers of h in the expansion (as we were 
expecting). Second, the constant term (i.e. the 0-th wave invariant), which corresponds to the term l = j = 
in the sum, equals 

a (t) = TrU (t) = TT(-g*_ )* \ , = TT 

j^wnukt' (2w fc tan^)* ^isin^ 

And third (using for j > 1 the coefficient of ft? in (|36|) equals 

k=l Sm 2 ;=i •'O ^0 JO 

The sum goes only up to 2j because if l > 2j then 2(1 + j) < 31 and Pj+j6;(0) = 0. 
This proves the first two parts of Theorem 11.11 

2.6. Calculations of the wave invariants and the proof of Theorem 11.11 3. In this section we try 
to calculate the wave invariants cij(t) from the formulas {7}. First of all, let us investigate how the terms 
with highest order of derivatives appear in dj(t). Because bi is the product of / copies of Wn functions, and 
because we have to put at least 3 derivatives on each Wn to obtain non-zero terms, the highest possible order 
of derivatives that can appear in Pj + ib[(0), is 2(j + l) — 3(1 — 1) = 2j — l + 3. This implies that, because 
in the sum {7} we have 1 < l < 2j, the highest order of derivatives in aj(t) is 2j + 2 and those derivatives 
are produced by the term corresponding to l = 1, i.e. Pj + ibi(0). The formula ([7]) also shows that aj(t) is a 
polynomial of degree 2j. The term with the highest polynomial order is the one with / = 2j, i.e. i"3j&2.j(0) 
(which has the lowest order of derivatives) and the term Pj + ibi(0) is the linear term of the polynomial. Now 
let us calculate Pj+i&i(x, z, £ ) and prove Theorem ll.11 3. 

By {37}, 

(38) P > +lbl - 2^.( ] + l)\ ^ A ^ ^ dr^.dr^ 

ri,...,r 2 j+2EAi J 



where here by ([20 



Also by 



&i = Wh ( C0SUJkS z k _ sinuj kS ^ k + ^ smuj k (t-s)+smuj k s k 
2 uj k sinw/jt 




H- 1 = 

V o I 

1 k k z K f K £ K Z K eK£K 

Hence the only non-zero entries of are the ones of the form hf x , h 1 = h\ , and h\ . Now we let 

i x k x k = the number of times hf x appears in h\ ir ' 2 ...h r 1 23+ir23+2 in (|38|) . 

i z t^fc = the number of times h\ ^ appears in h r 1 lT2 .../i^ 2j+ir23+2 in {38]) , 

i^*** = the number of times /i^ appears in /i^ ir2 .../i^ 2;,+ir23+2 in {38}. 
By applying these notations to {35}, {20} we get 
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i-U+D Q + l)!2^^ ^ -cot^ '^ 

• /i \ . • \ 2i t i' ■ / • \ % fc h t fr 



nj f ano^t - s) + sm^fcs \ * æ / cos uj k s y*H" ( -smuj k s \ * ' 4t I 2 j+2 w 
il sffi^ J J l^^J ^2 Q1 ,... 2Q „H/, 



X 

k=l 

where a k = i x k x k + *z fc f fc + > f° r k = 1, ...,n 



Ncxt wc write the above big sum as 



y (i + 1)' (A) _ y C? + l)! y Ty « fc ! 

So the coefficient of D 2 /^ 2 2a Wft in Pj+\b\, equals 



i-U+V (-iy +1 tt / 1 u kt (smu k (t- s) + sinw fc s N 2 



9 Cot_ ^- I : 7 ) - cos uj k s smuk s + cotw fc isin WfcS 



2-' +1 (n«fcO(IIWfc*) Al ^2 2 V sin Wfc i 
Now we observe that the term in the parenthcsis simplifies to 



1 cjfei / sin<x>fe(t — s) + sina->fcs\ . . 2 1 Wfct 
— cot — cosa>fcSsmwfcS + cotWfctsm w^s = — cot . 

2 2 V sinw/cf / 2 2 



So we get 



V ' |o|=j+l V 7 

Finally, by plugging (x, z, £ ) = into equation ([55)) and applying it to 0, we get ©. This fmishes the 
proof of Theorem ll.11 3. 

For future reference let us highlight the equation we just proved 

ri,...,r 2 3+ 2 e-4i 3+ |5|=i+l v 7 

where 

ly = ]y( C0StJfc5 ~fc _ sintJfcS ^ + ^ sintJfc(f-5)+sina;fcS fc 
2 cjfc sinwfet 

2.7. Calculations of J^ 1 Pj + 2&2(0), and the proof of Theorem 11.21 Throughout this section wc 
assume that V is of the form ©. Hence, the only non-zero Taylor coefficients are of the form D 2 g + V(0), 

or ft,^). where = (0, ...,0, 1). 

We notice that based on our discussion in the previous section, the Taylor coefficients of order 2j + l appear 
in Jq Jq 1 -P/+2&2(0), and they are of the form D 2 J +1 V(0)D~V(0). Therefore we look for the coefficients of 
the data 

(41) { D l^V(0)Dl s V(0y, \a\ = j - l}. 

in the expansion of aj(t). 
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PROPOSITION 2.4. In the expansion of a^t), the coefficient of the data f 2Æ+3e ^(0)^3e V(0), \a\ = j — l, 
is 



Therefore 



(2i)i+ 2 a\ \2u) 2) ^ l «„ + l A 2w„ 2' 9w 



(43) «i(*) = 7^4- Y U^^A S D% +2 V(Q) 

y J 3 ' (2i)i +1 ^ olV2u> 2' 2a w 



(2i)i 

\S\=j+l 



C2{n) 

' (2i)3+ 2 



{a polynomial of Taylor coefficients of order < 2j}. 



Before we prove Proposition [231 let us show how to use this proposition to prove Theorem ll.2l 

PROOF of Theorem 11.21 First of all, we prove that for all a, the functions 

(cot|i) Q , 

are linear ly independent over C. To show this we denne 

cot : (0,tt)" — > R", 
cot(æi, ...,x n ) = (cot(xi), cot(x„)). 

Because ujk are linearly independent over Q, the set {{^t, ^t) +irZ n ; ieM}fl (0, 7r) n is dense in (0, n) n . 
Since cot is a homeomorphism and is 7r-periodic, we conclude that the set {(cot(^i) , cot(^t); t E R} is 
dense in R™. Now assume 

E, LO \ ol 

cs{ cot —t) =0. 

Since {(cot(^t), cot(^t); t e R} is dense in R", we get 

y cs x 3 = o, 

c? 

for every X — {X\, ...,X n ) e R™. But the monomials X a are linearly independent over C. So cg = 0. 

Ncxt wc argue inductively to recover the Taylor coefficients of V from the wave invariants. Since 

n 1 

a ° (t) = n «s£f' 

we can recover Ilfc=i sm i 2 > anc ^ therefore we can recover {u>k} up to a permutation. This can be seen by 
Taylor expanding n&=i sm ii 2~- ^ e ^ x * ms P erm utation and we move on to recover the third order Taylor 
coefficient Z??- V(Q). This term appears first in a\(t). By Proposition l2.41 we have 



I a I — 2 

+ «<">I^(R ( ^ co '"r )J + R)( D| »^ (0) )'' 

+{a rational function of LOk}- 
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Now since the functions {(cot ^t) a }\s\=2 an d (3^(2^7" co * i %^) 2 + g^r) are linearly independent over C, 
we can therefore recover the data {-D 2( j^(0)}|5| =2 an d {D$g V^O) 2 } from ai(t). So we have determined 
the third order term D^g V(0) up to a minus sign from the first invariant ai(t). This choice of minus sign 
corresponds to a reflection. We fix this reflection and we move on to determine the higher order Taylor 
coefhcients inductively. 
Next we assume 

D le n V (°) and that we know a11 the Taylor coefhcients DJV(0) with m < 2j. We 

wish to determine the data {-E^a+Se V(0)}\s\=j-i an d {D^^V (0)}\s\=j+i^ from the wave invariant dj(t). 
At this point we use Proposition [23 and to hnish the proof of Theorem ll.ll we have to show that the set of 
functions 

{(„,§,)*; W -j + 1 }u{M|r- (^(^K^ + 4); l-l-J-1}. 

are linearly independent over C. But this is clear from our discussion at the beginning of proof. 

Proof of Proposition ^. 41 As we mentioned at the beginning of Section 2.7, the data -D^a+lie V(0)D\ s V(0), 
\a\ = j — l, appears first in o, (i) and it is a part of the term J Q f Q Sl Pj +2 b 2 (0). So let us calculate those terms 
in the expansion of Pj +2 b 2 (0) which contain D^^^g V(0)D^ s V(0). By ([20]) . since here l = 2, we have 

h(si,S2,x,x,zi,Z2,£i,&) = WiW 2 , where, 
(44) Wx = W( cos ^ Sl (gj 8 + 4) - sjR ^Ci + ( Si "^( t -^J+ Sin ^ s i ) x fc) i 

w _ w/-coso)fcS2 z fc _ sin^ fc s 2 cfc , / sin io k (t-s 2 )+sin u k s 2 \ fc\ 



Also from (l34l) we have 



(45) 



V 





















-I 


-I 











-I 


-I 





— D(uj cot(a;i)) 


—D(uj cot(uJi)) 


-I 


-I 


—D(uj cot(cJi)) 


— D(uj cot(aJi)) 



/ 5nx5n 



By fl37]) and (gH), P J+2 b 2 {0) = ^77 #2, where 5 2 is the following sum 



(46) 



S 2 = h r 2 ^...h r 2 2 ^ +i (W 1 W 2 ) ri .. 

ri,---,r 2j+i £A 2 



r 2 j+4 



(0), 



where A 2 = {x k , zj% z 2 , £1 , £ 2 and f° r every r, r' G -4 2 , ^ 2 r ^ s * ne (r, r')-entry of the matrix in (|45[) . 
We would like to separate out those terms in S 2 which include D^^g V (0) D\g V (0) . To do this, from the 
total number 2j + 4 derivatives that we want to apply to WiW 2 , we have to put 3 of them on W\ (or W 2 ) 
and put 2j + 1 of them on W 2 (or W\ respectively) . These combinations fit into one of the following two 
different forms 



(47) 



si = J2 h 2 

ri,...,r 2j+4 eA 2 



+ 3l*2i + 4 



(Wi) rir2r3 (W 2 ) 



r4,r$...r 2 j+4 



(0). 
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There are 2{j + + 2) terms of this form in the expansion of SV 

(48) S 2 2 = ^ ir2 -^2 2j+3r2j+4 (^ 1 ) rir3l . 5 (^ 2 ) r2 , r4 , r6 , r7 ... I . 2j+4 (0). 

There are 2 3 ( ^ terms of this form in the expansion of S 2 ■ 

Now, we calculate the sums S 2 and Sf. 

2.8. Calculation of S 2 . We rewrite S 2 as 

så= E K^hr^ E ^ 5r6 -^ 2j+3r23+4 (^2) r5 ..., 2J+4 ) r4 (^ 1 ) nr2r3 (o). 

ri,...,r 4 r 5 ,...,r 2 j+4 

Then from the definition of W2 in (|44| and also from (|45|) it is clear that we can apply (j40|) to the sum in 
the big parenthesis above. Hence we get 

(49) si = u + iy. E h.{^ coi ^) a ( ^ ^ ir2 ^ 3r4 (^ a w 2 ) r4 (Wi) rir2r3 )(o). 

|S|=j ' ri,...,r 4 

This reduces the calculation of S\ to calculating the small sum 

a\= E ^^(^Wwu^^o), (^ 2 = i? 2 2 ^ 2 ). 

ri ,. . . ,r4 

Computation of the sum A\ is straight forward and we omit writing the details of this computation. Using 
Maple, we obtain 

If we plug this into (|49|) . after a change of variable a„ — ► a n + 1 in indices, we get 

(so, f E ^(g cot |)^_^ )( _± cot ^ )2Da+k v,o )Dk n».. 

|a|=J-l 

2.9. Calculation of Sf . We rewrite S"! as 

sf= e ^ ir2 ^ 3r4 ^ 2 5r6 ( E ^ 7r8 -^ 23+3r2j+4 (^ 2 ), 7 ... I . 2j+4 ) (mw 5 (o). 

ri,...,r 6 r 7 ,...,r 2 j+4 

Again from it is clear that we can apply (I40[) to the sum in the big parenthesis above. So 

(si) sf = (j+i)! e é(2§ cot y) a ( £ wtw^oyiu^o). 

a|=j-l ' ri,...,r 6 

So we need to compute 

A%= E ^ ir2 ^ 3r4 ^ 5r6 (W2) r2 ,r 4) re(^l)r ir3 r 5 (0), (W2 = D^ 2 W 2 )- 

ri,...,r 6 

Using Maple 

/ f 1 ^ æt ^) 2 " Y^l) ^k^r) (0). 

If we plug this into (fSTj) we get 
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(52) 

fj" **,*-(,+!)! Y. ^(^»4)°(-4 < ^ æt ^> 2 -T^> D -- l ' (0)D - K<0) - 

|S|=i— i 

We note that the part of the expansion of J Q f Q Sl Pj+2&2(0) which contains the data -D^S+Se ^(0) -^le 
equals 

£L {2l] + w + 1) j;£ « + *('J s >jff si). 

Finally, by applying equations ([50)1 and (152j) to this we obtain 

3. Appendix A 

In this appendix we prove Lemma 12.11 



Proof. First of all we would like to change the function 9 slightly by rescaling it. We choose < r < 2e 
so that h 1 ~ T = o(fr}~ 2e ). Then we defme 

(53) Q h (x) := e(^). 

Thus Ø/j G Cq°([0, oo)) is supported in the interval 1% = [0, ft T <J]. In Appendix B, using min-max principle 
we show that 

Tr{Q(H)e^ É ) = Tr{Q h {H)e^ É ) + 0(ft°°) = Tr{e=£ É ) + 0(h°°). 
Hence to prove the lemma it is enough to show 

Tr{Q(P)e^ Lp ) = Tr(Q h {H)e = ^ A ) + 0(ff°). 
To prove this identity we use the WKB construction of the kernel of the operators <d(P)e~ p and 
Qf l {H)e~n~ H and make a compression between them. 

3.1. WKB construction for Q{P)e~^ p . In |DSj| , chapter 10, a WKB construction is made for ®(P)e~~^ p 
for symbols P in the symbol class <5q(1) which are independent of H or of the form P(x,t;,h) ~ Pq(x,£ i ) + 
hPi(x, £) + where Pj S S® are independent of h (but not for symbols H = H(x, £, Ti) S S$ o ). 

It is shown that we can approximate @(P)e~i>~ p for small time t, say i G (— fo^o), by a fourier integral 
operator of the form 

U P {t)u{x) = {2nh)- n [ [ e l(vp(Lx - 7l) - y - r ' )/h bp(t,x,y,r 1 ,h)u(y)dydr 1 , 



where bp € C°°((— to, to)', S(l)) have uniformly compact support in (x,y,rj), and ipp is real, smooth and is 
defined near the support of bp. The functions tpp and bp are found in such a way that for all t € (— to, to) 

\\e(P)e^ p -u P (t)\\ tr = o(h°°). 

Let us briefly review this construction, made in |DSj| . First of all, in Chapter 8, Theorem 8.7, it is proved 
that for every symbol P £ Sq(1), we have 0(P) = Op w (ap(x,£,,h)) for some ap(x, £, H) G Sq(1), where here 
P and Op w (ap(x,£ l ,H)) are respectively the Weyl quantization of P and ap(x,£,H). It is also shown that 
ap ~ ap j o(a;,C) + hap^x, £) + ... for some ap t j(x, £) S «5q(1). The idea of proof is as follows. In Theorem 
8.1 of |DSj| it is shown that if € Cq°(R), and if € Cq(C) is an almost analytic extension of (i.e. 
8Q(z) = 0(|9tø|°°)), then 



(54) 
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Then it is verified that for some symbol r(x, £, z; fe), we have {z — P) 1 = Op w (r(x, £, z; fe)). By symbolic 
calculus, one can find a formal asymptotic expansion of the form 

r ( x £ z - fe) ~ -J— + h qi[x > g ' Z) + fe 2 ^ ^ Z) + 

by formally solving Op w (r(x, £, z; fe))j)a(,2 — P) = (z — P)$nOp u '{r(x, £, z; fe)) = 1. We can see that qj(x, £, z) 
are polynomials in z with smooth coefficients. Finally it is shown that 0(P) = Op w (ap(x 1 ^H)), where 
op £ is given by 

ap(x, £, fe) = — / 90(^)r(x, £, h)L(dz). 
By the above asymptotic expansion for r(x, £, z; fe) one obtains an asymptotic ap ~ ap.o + feop,i + where 



(55) ^ ^ ^ ljé(z) { q j[ X ^ 2 Z J +1 L(dz) = J_d^( qj (x,ttMt))\t=P {x , v )- 

Then, again in Chapter 10 of |DSj| , it is shown that <pp(t, x, rj) and bp(t, x, y, ry, fe) satisfy 

(56) dt<pp(t,x, rf) +P(x, d x (pp(t,x,r))) = 0, ipp\t=o = x.t], 
bp~bp,o + hb P:1 + ... 1 bp >j = b Ptj {t,x,y,ri) GC°°((-to,to); 

where 

{9*6p,j + (d x (p P , d x b P<j ) + \A x ip P . bp.-j = -±A x b PJ -i, j > 0, (bp,-i = 0), 
bp,j\t=o = ip{x,'q)ap^{^-,T])il){y,r l ). 

In (|57[) . ap.j is given by (|55[) and ^(x,^) is any function which equals 1 in a neighborhood of P _1 (7) 
where / = [0, S] is, as before, the range of our low-lying eigenvalues and where is supported. 

There exists a similar construction for Qj l (H)e = ^ H , except here H G Sg o . 

3.2. WKB construction for Q h (H)e^ É . Since in (JT3J), i? = H(x,£,K) G S# D , with <5 = § - e, we can 
not simply use the construction in |DSj| mentioned above. Here in two lemmas we show that the same 
construction works for the operator Qn(H)e~ H . We will closely follow the proofs in |DSj| . 

Lemma 3.1. 

1) Let Qfr be given by h53\) and H G S® o by U3\) . Then for some an G S$ we have Q(H) = 
Op w {auix, £, fe)). Moreover au(x, £, fe) ~ a^ot 1 , £■ fe) + faa>H,i{%, £, fe) + where aii,j(x, £, fe) G Sg o 
is given by 



(58) a HJ = ^J q -jøriør L ( dz ) = j^d?( qH Ax,tt,h)e h (t))\t=H {x ^H)- 

2) Choose c such that < c < min{l, a;|}^ =1 < max{l, w 2 }J? =1 < K Let ^)%{x,rj) be a function 
in Cg°(R 2n ) n S^ o (R 2n ) which is supported in the ball {x 2 + if < 4c _1 fe 1_T <5} and equals 1 in a 
neighborhood o/-ff _1 (i), where In — [0, fr 1 ~ T å] (In is where &n is supported). Then 

(59) e h (H)u(x) = (2irh)~ n J J e^^^^M*, v)a H (^-,V, VMv, vHy)dyd V + K{K)u(x), 
where \\K{h)\\ tr = 0(fe°°). 
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Proof OF Lemma 13.11 Since H G Sg and 5q = \ — e < |, the symbolic calculus mentioned in the last 
section can be followed similarly to prove Lemma l3"7Tl l. It is also easy to check that in ([IS]), a H ,j G S$ o . The 
second part of the Lemma is stated in |DSj| , equation 10.1, for the case P G Sg. The same argument works 
for H G Sg , precisely because the factor h N on the right hand side of the inequality in Proposition 9.5 of 
|DSj| changes to h N ~ s ° a . Thus the discussion on pages 115 — 116 still follows. 

Lemma 3.2. For every t in some small interval (—to, to), there exist functions (fn(t, x, r\, K) and bu(t, x, y, rj, K) 
such that the operator Un(t) defined by 

(60) U H (t)u(x) = (27T/1)-" j J e iW H {t,x^n)-y. v )/h bH ^ X) yi ^ h ) u ( y ) dyd ^ 

satisfi.es 

\\Q H (H)e^ A -U H (t)\\ tr = 0(h°°). 



Moreover, we can choose ipn and bu such that 

1) ifH satisfies the eikonal equation 

(61) d f (p H (t, x,n, K) + H(x,d x ip H (t,x,r), K)) = 0, ip H \ t=0 = x.n. 

This equation can be solved in (—to, to) x {x 2 + rj 2 < Ch 1 ~ T 6} where C is an arbitrary constant. In 
faet ipu is independent of h in this domain. (Only the domain ofcfH depends on h. See \64\) -) 

2) For all t G (—to, to), we have bu(t, x,y,n,h) G Sg o with supp bu C {x 2 + n 2 , y 2 + rj 2 < C\h 1 ~ T 5} for 
some constant C\. Also bu has an asymptotic expansion of the form 

(62) b H ~b H ,a+f&H,i + -.., b H , 3 ; = b H ,j(t,x,y,r),h) G C°°((-t ,to); S° So (l)), 
and the functions bnj satisfy the transport equations 

{dtbnj + (d x <fiH,dxbH.j) + ^A x tp H .b H ,j = -^A x b H ,j-i, j > 0, (b H .-i = 0), 
x 
bH, 3 \t=a = ^n(x,v) a H,j( S: r-,'n>K)'*Ph(y,'n), 

where in H63\) we let ipfr(x,r]) be a function in C x (M. 2n ) n Sj^R 2 ™) which is supported in the ball 
{x 2 + n 2 < 4ch 1 ~ T å} and equals 1 in a neighborhood of H^ 1 ^^), where Ih — [0, h 1 ~ T 6]. Here c is 
defined in Lemma \3.1\ 2. Also in \6S\) , the functions an j are defined by L58\) . 

4) 

(64) ip H (t,x,n,h) = ip P (t,x,if) on {x 2 + n 2 , y 2 + i] 2 < Cih 1 ~ T 5} D supp (b H (x, y, n, H)). 

5) 

(65) b H j(t,x,y,r),H) = b PJ (t,x,y,n) on {x 2 + rj 2 , y 2 + n 2 < ch 1 ~ T 5}. 

Proof of Lemma 13.21 First of all we assume Un(t) is given by (|60[) and we try to solve the equation 

\\(!fd t +H)U H (t)\\ tr = 0(h™), 
U H (0) = &h(H). 

for ipH and bn, for small time t. 
Using (|59|) . this leads us to 
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e -*>PH/K^ dt + H)(e^"/ h b H ) G C°°((-t ,to); ^(l)), 

b\ t =o = ipn(x, 77)0^(^,77, h)ipH(y, v)- 
We choose the phase function = yjfj(i, x, 77, h) to satisfy the eikonal equation (f6Tj) . This equation can 
be solved in a neighborhood of the support of bjj, for small time t G (— to, to) with to independent of Ti. 
Let us explain how to solve this equation. We let (x(t,z,rj;h),^(t,z,r];h)) be the solution to the Hamilton 
equation 

)d t x — d$H(x,(;, h) — £, x(0, z, 77; h) = z 

■ 
dtt = -d x H(x, t K) = -d x V H (x), £(0, z, 77; h) = r, 

We can show that (see section 4 of |Chj ) there exists to independent of h such that for all \t\ < to we have 

!\d z x(t,z,r);K) - I\ < \, \d v x(t,z, r);K)\ < \ 
■ 
\d z £(t,z,rr,h)\ < \, \d n £{t,z,r);h)-I\ < \ 

We can choose to independent of H, precisely because in equation 4.4 of |Chj we have a uniform bound in h 
for Hess(V?i(x)). Now, we define 

A : {z,rj) 1 — > (x(t,z,rj;h),ri). 

It is easy to see that A(0,0) = (0,0). This is because if (z,T)) = (0,0) then H(x,£) = H(z,7]) = 0. By ([13]) 
and ©, and W{x) = 0(|x| 3 ), we can see that ff(x,£) = implies (x(t, 0, 0; H),£(t, 0, 0; Ti)) = (0,0). On the 
other hand from (|67[) we have | < \d z x(t, z, 77; h)\ < |. Therefore Å is invertible in a neighborhood of origin. 
We define the inverse function by 

X^ 1 (x,t]) = (z(t,x,rj; H),rj), 
which is defined in a neighborhood of (x, 77) = (0, 0). Then we have 

//* 1 

(68) <pit(t, x,rj, h) = z(t,x,Tj; Ti). 77 + J ^{s, z{t, x, 77; Ti), 77; K)\ 2 - V h (x(s, z(t, x, 77; K), 77; H))ds, 

A similar formula holds for ipp except in (l66|) H should be replaced by P and in (|68|) V% by V. It is known 
that the eikonal equation for ipp can be solved near supp&p, for small time t G (— to, to) (Of course to is 
independent of h). Now, we want to show that 

(69) ip H (t,x,r),ti) = <p P (t,x,T)) in (-t , t ) x {x 2 + r? < Ch 1 ~ T 6}. 

Let (x, 77) be in {x 2 +r/ 2 < Ch 1 ~ T 5}. First, we show that \z(t, x, 77; h)\ < SC^h^S^ . Because z(t, 0, 0; Ti) = 0, 
by Fundamental Theorem of Calculus we have 

\z(t,x,ri;h)\ < (\x\ + |r7|)sup{(|9 æ | + \d v \){z{t, x, 77; h))}. 

From x(t, z(t, x, ri; K), 77; Ti) = x, we get 

d v z — —(d z x)~~ d v x. 

Thus by ((6Tj> , \d x z\ + \d v z\ < 4. Hence \z(t, x, 77; K)\ < 4(|x| + |?7|) < SC^H^å^. This implies that for 
all \t\ < t , {x(s, z(t, x, 77; h), 77; Ti), £(s, z(t, x, 77; Ti), 77; Ti)) will stay in a ball of radius 0(h 1 ~ T ) centered at the 
origin. On the other hand, by definition (fT51) . P and H agree in the ball {x 2 + rj 2 < jh 1-25 } and r < 2e. So 
for all i, s G (— 1 , t ) and (x, 77) G {x 2 + 7/ 2 < Ch 1 ~ T å} we have 



z P {t,x,rj) = Z(t,x,r);H), 
(70) x P (s, z P (f , x, 77), 7/) = x(s, z(f , x, 77; Ti), 77; Ti), 

£p(s, z p(t, x > v),v) = z (t, x, 77; Ti), 77; Ti), 
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where Zp(t,x,rf), Xp(s,Zp(t,x,rf),rj) and £_p(s, zp(t, x,rj),r]) are corresponded to the Hamilton flow of P. 
Hence by and a similar formula for ipp, we have (joU)) . This also shows that we can solve (|5Tj) in 
(-io ! io)x{x 2 +r; 2 <Cft 1 - r 5}. 

To find bf{ we assume it is of the form (|62|) and we search for functions bn.j such that e~ lVH ^ h (jdt + 
H)(é lipH / h bH) ~ 0. After some straightforward calculations and using the eikonal equation for we obtain 
the so called transport equations ([63]) . Now let us solve the transport equations inductively (see [Ch] ) . 

In |Chj it is shown that the solutions to the transport equation ([63]) are given by 

b H . (t, x, y, ry, K) = J~?(t, x, ry, H)b H . (Q, z(t, x, ry; K), ry; h),y, ry, Ti) 

( 71 ) b H .j (t, x, y, ry, h) = J~? {t, x, ry, %) (b H ,j (0, z(i, x, 77; Ti), ry; Ti), y, ry, Ti) 

^ 5 Jo J ^ ( s ' æ ' ^ h ) ^ b H,j-l (s, æ(s, z(i, x, 77; ft), ry; Ti), y, ry, Ti)ds) . 

where 

J(t, x, ry, Ti) = det(9 x #(t, x, ry; Ti)) -1 . 
Now, we notice by the assumption on we have supp(6jj,j(0, x, y, ry; Ti)) C {x 2 + ry 2 , x 2 +ry 2 < 4c _1 Ti 1 ~ r (5}. 
So by our previous discussion on z(i, æ, ry, Ti), we can argue inductively that for all t G (— <o, io)j supp(bfj/,j) C 
{x 2 + rj 2 7 y 2 + ry 2 < Cifr 1-1 ^} for some constant C\. Since bH,j\t-o é Sj , we can also see inductively from 
(J7TJ) that S S 1 ^. Finally, Borel's theorem produces a compactly supported amplitude bn € 5° from the 
compactly supported functions bu j € S g . This finishes the proof of items 1 — 3 in Lemma 13.21 

Now we give proofs for items 4, 5 in Lemma 13.21 

By choosing C > C\, equation is clearly true from ([BU]) . Next we prove that equation ([US]) holds. 
Using (|55|) and ([55]) . and because P and Tf agree in the ball {x 2 +ry 2 < ifi, 1 ~ 2e }, we observe that the functions 
cip,j(x, ry) and ajj,j(x, £, Ti) agree in this ball. Therefore, because suppen (x, ry) C {x 2 + rj 2 < 4c~ 1 Ti 1 ~ r (5} and 
ijj h = 1 in {x 2 + rj 2 < ch 1 - T 5}, by ((571) and |B3J 

b Htj (0,x,y,ri,h) = b Pij (0, æ, y, ?y) on {(2, ry, ry); æ 2 +?y 2 , y 2 + rj 2 < ch 1 ~ T 5). 

This proves ([65]) only at /j = 0. But by applying ([70]) to ([71]) and a similar formula for bp, we get (|65|). This 
finishes the proof of Lemma |3~21 

To finish the proof of Lemma l2~T] we have to show that for t sumciently small TrUn(t) = TrUp(t)+0(h°°), 
or equivalently 



e i( VS {t,x,r,,H)-x.r,)/h bH ^ ^ ^ ^ ty dxdrj _ / / e i( w (t,i fl )-x.,)/fi 6p ^ ^ ^ ^ fi)d:rd?y + 0(ft°°). 

By (|64p . the phase function (y9# of the double integral on the left hand side equals ipp on the support of the 
amplitude bn, so ipn is independent of h in this domain. Now, if t £ (0,in) where to is smaller than the 
smallest non-zero period of the flows of P and H respectively in the energy balls 

{(x,r,)\ H{x,fj) < 5h l - T } c {(x,r,)\ P(x,r,) < S}, 

then for every such t, (x, 77) = (0, 0) is the only critical point of the phase functions (pH(t, x, ry, K) — x.rj and 
<pp(t, x, ry) — x.t) in these energy balls. 

Obviously both integrals in the equation above are convergent because their amplitudes are compactly 
supported. But the question is whether or not we can apply the stationary phase lemma to these integrals 
around their unique non-degenerate critical points. By Lemma 13.21 the phase functions ipu and tpp are 
independent of h on the support of their corresponding amplitudes. Hence ipH,<fip E S® on suppfo# and 
supp bp respectively. On the other hand bn(t, x, x, 77, h) E S® , Sq < h; and bp(t, x, x, rj, K) E Sq . These facts 
can be used to get the required estimates for the remainder term in the stationary phase lemma (for an 
estimate for the remainder term of the stationary phase lemma, see for example Proposition 5.2 of |DSj| ). 




INVERSE SPECTRAL PROBLEMS FOR SCHRODINGER OPERATORS 



21 



Finally, by (|64|) and ([55]) it is obvious that the integrals above must have the same stationary phase 
expansions. 

4. Appendix B 

In this appendix we prove Lemma \2. 21 In faet we prove that if Qn is given by (l53|) then in the sense of 
tempered distributions 

(72) Tr{Q h (H)e^ É ) = Tr{e^ k ) + 0(h°°). 

Proof of Lemma 12.21 follows similarly. 

We will use the min-max principle. 

Min-max principle. Let H be a self-adjoint operator that is bounded from below, i.e. H > cl, with 
purely discrete Spectrum {Ej}°°^ a . Then 



(73) 



Ej = sup 



inf 



ip u ..., Vn -! ^ g D(H); \\tp\\ = 1 

tp G span(<y9i, ifin-i) 1 - 



As before we put H = -±h 2 A+V h (x) = -±h 2 A+± Y2=i 0J 2 k x 2 k +W h {x) , and H = -Aft 2 A+± ^Li 



Then if we let C = \\W h (x) 



i \ , we have 

i L-(R"x(0,/io)J ' 

(V>, #oV0 — C < (i/j, Htp) < {xjj, H ^) + C, 



and therefore by applying the min-max principle to the operators H and Hq we get 
(74) E° (fi) - C < E 3 (h) < E° {h)+C. 

Notice we have explicit formulas for the eigenvalues Ej(H) of Hq. They are given by the lattice points in the 
first quadrant of R™. More precisely 



n 

a(H ) = {Ei(h) = hJ^Mlk + -); lk G Z^ }. 



fc=i 



Since in the sense of tempered distributions 

Tr{Q n {H)e = ^ 6 ) = Tr( X[aM ^] (H)e^ å ) + 0{h°°)- (see for example [KD]), 
to prove (|72|). it is clearly enough to show that for every tp in S(M.) 



E ØF^p) = o(»~). 

{j;B J (fi)>5fti-T} 

Since ip is in 5(M), for every p > there exists a constant C p such that 



\p(x)\<C p \x\-P. 



Hence by (0 



< Cr 



E°Ah)-C 



Again using ([74]) and because C = ||Wk(a;)|| T /_ , n , ,\ < ^4^~ 3e < f ft 1 r wc get 

L°° ^JR n x (0,/iq)J 



< c„< 



åh 1 



c 



5h l - T - 2C" 



E°ÅK) 



<2<X 



Ej (fy 



for Ej(h) > Sh 1 
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Now let m be an arbitrary positive integer. So in order to prove the lemma it is enough to find a uniform 
bound for 

" 1 

By applying the geometric-arithmetic mean value inequality we get 

n 1 

A(7i) < Tl^TT™ £ | JJ Wfe(7fe + _)|-P 

{7; E^(7*+^)> " 1 7- C } fc=1 

n 11 

<n- p E{( s " E K^+^r) n (E^(*+2 rp )}- 

{7kG z - ; "fc(7fc+;j)> } 

We claim for p large enough there is a uniform bound for the sum on the right hand side of the above 
inequality. It is clear that if p > 2 then the series J2 7 , l w fe'(7fc' + \)\~ p lSi convergent. Also if for some 7fc 

we have u)k{lk + 3) > &h „^~ C ; then because C = 0(h^~ 3e ), for small enough we have (uJk{^k + é)) 1 T > 

(^)Vrl. Thus 

{7fe ez>o ; Wk(7fc+ i)>«iz^ } 7* 

So if we choose p > max { — , 2}, then the sum on the right hand side is convergent and therefore we have a 
uniform bound for the sum on the left hand side and hence for A(K). This finishes the proof of (|72p . 
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